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A simple and very flexible variational approach to the out-of-equilibrium quantum dynamics in 
strongly correlated electron systems is introduced through a time-dependent Gutzwiller wavefunc- 
tion. As an application, we study the simple case of a sudden change of the interaction in the 
fermionic Hubbard model and find at the mean field level an extremely rich behaviour. In particu- 
lar, a dynamical transition between small and large quantum quench regimes is found to occur at 
half-filling, in accordance with the analysis of Eckstein et a7., Phys. Rev. Lett. 103, 056403 (2009), 
obtained by dynamical mean field theory, that turns into a crossover at any finite doping. 
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Introduction. Triggered by the enormous advances in 
the physics of ultra-cold atomic gases, [l| time dependent 
non equilibrium phenomena in strongly interacting quan- 
tum systems have recently become of greatest interest. 
The possibility of artificially engineering many-particle 
quantum states with tunable interactions and almost per- 
fect isolation from the environment gives the chance of 
probing directly in the time domain the quantum dy- 
namics following an external perturbation y. While 
early experiments focus mainly on bosonic systems [3| 
or fermionic condensates, [^] the recent experimental re- 
alization of a fermionic Mott insulator [5[ opens the way 
to investigate out-of-equilibrium phenomena in electron 
systems too.[6i] From a theoretical perspective, these ex- 
periments raise several intruiguing questions touching 
quantum dynamics at its roots. Indeed, when driven 
out of equilibrium, interacting quantum systems can dis- 
play peculiar dynamical behaviours or even be trapped 
into metastable configurations [7| that differ completely 
from their equilibrium counterpart. The simplest way 
one could imagine to induce a non trivial dynamics is 
through a so called quantum quench. Here the system is 
firstly prepared in the ground-state of some given Hamil- 
tonian "Hi , and then suddenly let evolved under the action 
of a new hamiltonian "Hy. Recently, quantum quenches 
have been the subject of a vast literature focusing on inte- 
grable systems, [^ one dimensional models, [9|] or systems 
close to a quantum criticality 10] The interest on these 
class of non equilibrium problems relies both on the dy- 
namics itself. 111 as well as on the long-time properties 
where the issue of thermalization or its lack of is still 
highly debated. l^jllSl For what concerns strongly corre- 
lated electrons in more than one dimension, the subject 
is still largely unexplored and progresses have been done 
only very recently. The single band fermionic Hubbard 
model is likely the simplest lattice model of correlated 
electrons emboding the competition between metallic and 
insulating behaviour driven by a local Hubbard repulsion 



U . Its Hamiltonian reads 
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In two pioneering works. |14l.ll5l| the response of a Fermi 
sea to a sudden switch-on of the Hubbard U has been 
studied in infinite dimensions using respectively the flow- 
equation method and the Dynamical Mean Field Theory 
(DMFT). Results suggest the existence of two different 
regimes in the real-time dynamics depending on the fi- 
nal interaction strength Uf. At weak coupling, [14| the 
systems is trapped at long-times into a quasi-stationary 
regime where correlations are more effective than in equi- 
librium. This pre-thermalization phenomenon has been 
confirmed by DMFT results, [15] which further indicate 
a true dynamical transition above a critical C//c towards 
another regime with pronounced oscillations in the dy- 
namics of physical quantities. These intriguing results 
have been so far restricted to a quench starting from a 
non interacting system (Ui — 0) and, more importantly, 
limited to rather short accessible time scales and weak 
quenches, thus leaving open many important issues. 

A simple and flexible approach, although less rigor- 
ous than e.g. DMFT, is thus extremely desirable and 
this is actually the aim of the present work. Specifi- 
cally, here we propose a variational approach to the out- 
of-equilibrium dynamics of a correlated electron system 
based on a proper extension of the Gutzwiller wavefunc- 
tion. We apply this technique to study the interaction 
quench in the Hubbard model, where we find a rich be- 
haviour featuring a transition in the real-time dynam- 
ics at a critical quench line Ufc (Ui), in accordance with 
Ref. [l5| at Ui = 0. Remarkably, a finite doping com- 
pletely washes out this transition, leaving behind only a 
crossover from weak to strong coupling. 

Variational Approach to Quantum Dynamics. Many 
basic concepts in the theory of strongly correlated sys- 
tems, like e.g. the Brinkman-Rice scenario for the Mott 



transition, have been originated from calculations based 
on a very simple and physically transparent variational 
approach introduced in the 60th's by Gutzwiller [l6| . 
This approach has been so far applied only at equilibrium 
or at most in the linear response regime. [171] but it turns 
out to be so flexible to allow for full out-of-equilibrium 
calculations. 

For simplicity, we assume initially a many-body wave- 
function |*o), which, for times i > 0, is let evolve with 
a Hamiltonian Ji that includes sizable on-site interac- 
tions. In the spirit of the Gutzwiller approach, we make 
the following variational ansatz for the time-dependent 
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where |<&(i)) is a time-dependent Slater determinant. 
Vi{t) is a hermitean operator that acts on the Hilbert 
space of site i and controls the weights of the local elec- 
tronic configurations. Si{t) is also hermitean and we as- 
sume it depends on some variables </),;„ (i) such that 
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where Oia is any local hermitean operator. Since ([2]) 
is just a variational ansatz, it does not solve the full 
Schroedinger equation. Our proposal is to determine the 
variational parameters by requiring: (i) that the Heisen- 
berg equations of motion of the local operators Oia are 
satisfied when averaging over ([2]); (ii) that the average 
energy E = {'^{t)\T-L |^(t)) is, as it should be, conserved 
during the evolution. Since, by definition. 
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where the equivalence is our variational assumption. 
Within the Gutzwiller approximation, [18J which is ex- 
act in the limit of infinite coordination lattices, E = 
{^{t)\'H*{t)\^{t)), where ■H*(t) is a non-interacting 
Hamiltonian that depends on all time-dependent vari- 
ational parameters defining Vi and Si. In general, these 
parameters can be expressed in terms of (j^ia and Oia ■ If 
we impose that |$(i)) is the solution of the Schroedinger 
equation, namely that —idt \^{t)) = 'H*(t) |$(t)) and fur- 
thermore that 



dt 



dE 

dOia 



(4) 



conservation of energy follows automatically. Therefore, 
4>ia and Oia act like conjugate variables and the energy 
functional E as their effective Hamiltonian. 

As a simple application of the above variational 
scheme, we assume Ti. to be the Hubbard model ([T]) at 
half-filling with U{t < 0) = [/, > and U{t > 0) = 
Uf > Ui, and furthermore we limit our analysis to ho- 
mogeneous paramagnetic wavefunctions. In the limit of 
infinite coordination lattices, one can compute exactly 
average values on the variational wavefunction provided 
the following conditions are imposed [ISJ, ll9| 
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We assume Vi{t) = X],^i=o ^i,7i{t)'Pi,n, where Vi^n is the 
projector at site i onto configurations with n = 0, . . . , 2 
electrons and Si{t) — J2n=o 4>i.n{t)'Pi^n, which implies 
that (t>i^n [t) plays the role of the conjugate variable of the 
occupation probability Pi_„ — (^(i)|7'i_„|\l'(t)). From the 
constraints above it follows that Pi^ = P, 2 and Pii = 
1 — 2Pj 2- We define Pi^2 = (1 ~ cos6'j)/4 and set (t)ifi = 
<f>i2 = 4>i while <f>ii = 0. Using 9i and (pi as variational 
parameters, one finds the average energy [18J 

^ = ^Y.^l-^os6,{t)) (5) 
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H.c.\<^{t)). One 



recognizes in ([5]) the mean field energy of an Ising model 
in a transverse field 
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where {a^) — cos 9i and (trf) = sin^i cos^^. This 
connection can be established rigorously at the varia- 
tional level, [20j and agrees with the .^2-slave-spin the- 
ory recently introduced. [2l|, |22| Therefore, it is not 
surprising that the equations of motion that we obtain 
through Eqs. ^ and (g]) are just those of the Ising model 
{dtcrf) = —i{[a°;,'Hi\) within mean field. Therefore, un- 
der the above assumption of homogeneous and param- 
agnetic wavefunctions, a quantum quench in the half- 
filled Hubbard model is equivalent, within the Gutzwiller 
variational scheme, to a quench in a Ising model in the 
presence of a transverse field. In particular, if |<I'(t)) is 
taken to be the half-filled Fermi sea, then Wij{t) = w 
and © is the conventional ferromagnetic Ising model 
with constant and uniform exchange w and transverse 
field C///4. Quantum quenches of the transverse field 
have been recently investigated in one-dimension J8| and 
on a fully connected lattice [23|. In the following we 
assume that the system is prepared in the metallic varia- 
tional wavefunction that optimizes ([1} with U = Ui < Uc, 
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FIG. 1: Left Panels: Gutzwiller mean field dynamics at iialf- 
filling for quasiparticle weight Zit) (black line) and double 
occupation D(i) (dashed red line) for quantum quenches from 
Ui = 0.25 to uj — 0.35 (top panel) and Uf = 1.25 (bottom 
panel). Right Panel: Period of oscillations at half filling and 
for a finite doping S j^ 0. Note the logarithmic singularity at 
ufc for (5 = (see main text) which is cut-off by finite doping. 



FIG. 2: Average double occupation D (top) and quasiparticle 
weight Z (bottom) as a function of uf at fixed Ui = 0.0, 0.5. 
We show results at half filling (full lines) that display a sharp 
transition at Ufc, as well as at finite doping (dashed lines) 
where only a crossover remains. We also plot the zero tem- 
perature equilibrium results for 5 = (red points). 



where Uc is the variational estimate of the Mott transi- 
tion. This corresponds to initial values 4>i(0) = and 
cos9i{0) = Ui/Uc = Ui for the coupled equations: 2(j) = 
Uc cos 6 cos^ 4>~Uf and 20 = Uc sin sin (p cos (/). We note 
that, apart from the trivial case in which Uf = Ui, these 
equations admit a non-trivial stationary solution 9 = 
and cos^ 4> = Uf/Uc = Uf, which is compatible with the 
initial conditions only when Uf = Ufc = (1 -I- Ui)/2. It 
turns out that Ufc identifies a dynamical critical point 
that separates two different regimes similarly to a simple 
pendulum. When Uf < Ufc, 2(f)(t) oscillates around the 
origin, while, for Uf > Ufc, it performs a cyclic motion 
around the whole circle. 

In order to characterize the different regimes, we 
focus on three physical quantities, the double occu- 
pancy D{t) = (1 — cos6'(i))/4, the quasiparticle residue 
Z{t) = sin^ 9{t) cos^ 4'{t) and their period of oscillation, 
T. While detailed calculations will be presented else- 
where [20|, in the rest of the paper we just sketch the 
results of the mean field dynamics. Let us start from 
the weak coupling side Ui < Uf < Ufc, see top panel 
in Fig. [U where both D{t) and Z{t) display small oscil- 
lations. Their amplitude and period increase with the 
strength of the quench Su = Uf — Ui, the latter reading 
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Uf > Ufc, a novel strong-coupling dynamical behaviour 
emerges. Here oscillations become faster, their period 
T = 4/C (1/fc) / yJufSu now decreases as a function of 5u. 
In particular, for Uf ^ Ui, we get T — — smoothly 
matching the atomic limit result. The oscillation ampli- 
tude of -D(i) decreases with Uf, which results into a frozen 
dynamics in the infinite quench limit, [7] while quasipar- 
ticle weight Z{t) still shows large oscillations even for 
M/ — >■ cx), mainly reflecting the unbounded dynamics of 
the phase (f){t). 

Remarkably, the weak and the strong coupling regimes 
are separated by a critical quench line Ufc at which mean- 
field dynamics exhibits exponential relaxation. Indeed, 
upon approaching this line from both sides, the period 

r diverges logarithmically, T ~ -/= log ( |„^.i„c| j ■ 

Right at criticality, Uf = Ufc, the mean field dynam- 
ics can be integrated exactly. The result gives D{t) = 
£'(0) (l — tgh^ (i/r*)); the double occupation relaxes ex- 
ponentially to D = pushing also Z — ?► 0, with a char- 
acteristic time scale r* = 2/y^Z{0) that increases upon 
approaching the Mott Insulator. 

We now turn to discuss long-time average properties 
of the Gutzwiller mean-field dynamics that we define 
as O = limt^ooj/n dt' 0{t'). The analytical expres- 
sions [20| of D and Z are shown in Fig. [21 At weak 



where /C (fc) is the complete elliptic integral coupling, we find D = D{0) [l + ? ( 



of the first kind with argument k^ 
(5u — > we find a linear increase T 

withTB =^TT/y/Z{0). 
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, where 

E (k) is the complete elliptic integral of the second kind 
and fc^ is the same as before. In addition, due to energy 
conservation, the knowledge of D completely fixes the av- 



Conversely, when quenching above the critical value, erage quasiparticle weight Z = Z{0) + 8uf {D — Z?(0)). 



It is interesting to consider first the small quench limit 
5u — >■ 0. We find that, given D{0) = (1 — Ui)/A the initial 
equilibrium value, D ~ D{0) — 6u/4 = (1 — u^)/4, namely 
tends to the equilibrium value corresponding to the final 
interaction. Hence the prethermalization result [IJ] for 
the quasiparticle weight Z immediately follows. Indeed, 
quenching from a non-interacting Fermi sea, Ui — hence 
Z(0) = 1, we find that the non-equilibrium Z is reduced 
twice more than its equilibrium value at Uf. 

For large quenches, the average double occupancy D 
increases as a function of the final interaction Uf,D — 
D(0) -\{uf- Ui) {l-E (1/fc) //C (l/fc)), eventually ap- 
proaches its initial value -D(O) for Uf ^ oo. A similar be- 
haviour is found for Z which is however further reduced 
by a factor 1/2 with respect to the initial value Z{0) due 
to the freely oscillating behaviour of the phase. We find 
therefore that, for large quantum quenches, the dynamics 
retains memory of the initial condition and thermaliza- 
tion is prevented by a dynamical blocking. 

Finally, for quenches close to u/c, both D and Z are 
very small, vanishing as l/logju/ — u/d on approach- 
ing the critical point. Namely, ujc not only signals a 
transition in the dynamics but also identifies the critical 
interaction at which the quenched system shows genuine 
Mott insulating behaviour. 

Away from half-filling, the dynamical equations be- 
come more cumbersome [20|. However, key features can 
be easily derived even without resorting to a numerical 
integration. In particular, we find that any finite dop- 
ing turns the half-filled dynamical critical point into a 
crossover. For instance, the logarithmic singularity of 
the oscillatory period T is cut-off by any finite doping, as 
shown in figure [1] As a consequence, the singular behav- 
ior of the average values across the half-filling transition 
is smoothed into a crossover at finite doping, see figure 

m 

Discussion. It is worth discussing the above results in 
light of those recently obtained by DMFT.^] Remark- 
ably, our variational ansatz ([2]) seems to catch many non 
trivial effects observed in DMFT. In particular the exis- 
tence of two different regimes separated by a real dynam- 
ical transition at Ufc, already suggested in [15|, clearly 
emerges from our mean field theory. 

We note however that the suppression of quantum fluc- 
tuations, which is at the ground of our results, give rise 
to an oversimplified periodical dynamics that lacks re- 
laxation. In this respect DMFT, which can treat exactly 
all local quantum fluctuations, works much better and 
turns these oscillations into a true relaxation. We guess 
that a similar result could be obtained from our varia- 
tional treatment, for example, by allowing fluctuations 
in the Fermi sea and treating the Quantum Ising Model 
beyond the simplest mean field level 2]|. Finally we 



role of the neglected quantum fluctuations is to provide 
dissippation for the two-level system described by the lo- 
cal Ising variable, the two dynamical regimes found by 
DMFT resemble those in the phase diagram of the spin- 
boson model. 24 1 

Conclusion. We introduced a variational approach to 
the out-of-equilibrium dynamics in strongly correlated 
electron systems. Using a time dependent Gutzwiller 
ansatz we address the problem of an interaction quench 
in the Hubbard model finding a dynamical transition at 
half-filling. Our results provide a simple and intuitive 
mean field theory for the quench dynamics in interacting 
Fermi systems. 
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notice that the Ising analogy provides a simple interpre- 
tation of the dynamical transition, at least from a local 
DMFT-like point of view. Indeed, if we assume that the 
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